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(57) ABSTRACT 

A pseudonise (PN) code acquisition scheme employs all 
elements of a smart antenna array and an adaptive threshold. 
The basic structure is the combination of a conventional PN 
correlation searcher, an adaptive beamformer, and an adap- 
tive threshold setting circuit. During each observation 
interval, the adaptive beamformer adaptively updates the 
weight vector for the smart antenna elements using the 
accumulated received signal despread with trial PN code 
phase error as input, preferably at the chip rate. A spatially 
correlated signal is then formed by weighting an accumu- 
lated value of the signal received by each antenna in the 
array over the observation period with the corresponding 
final weight of smart antenna weight vector as calculated by 
the adaptive beamformer. This spatially correlated signal is 
then compared to a threshold to determine whether PN code 
acquisition has occurred. In preferred embodiments, an 
adaptive threshold setting algorithm is employed to adapt to 
varying environment for efficient PN code acquisition. The 
adaptive threshold setting circuit accumulates the signal 
received by each element over the observation period, 
multiplies the accumulated values by the updated weight 
vector and combines the weighted values to an average 
estimated power. The average estimated power is then 
employed to scale a fixed reference threshold to create the 
adaptive threshold. 
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PN CODE ACQUISITION WITH ADAPTIVE 
ANTENNA ARRAY AND ADAPTIVE 
THRESHOLD FOR CDMA WIRELESS 
COMMUNICATIONS 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to wireless telecommunica- 
tions generally, and more particularly to the application of a 
smart antenna to PN code acquisition in a code division 
multiple access (CDMA) wireless communications system. 

2. Description of the Related Art 

In third generation (3G) wireless communications 
systems, particularly wide band (W)-CDMA systems such 
as those described in the 3rd Generation Partnership Project, 
"Physical Channels and Mapping of Transport Channels 
onto Physical Channels (FDD)/' 3GPP Technical 
Specification, TS25.211, v3.2.0, March 2000; 3rd Genera- 
tion Partnership Project, "Spreading and Modulation 
(FDD)," 3GPP Technical Specification, TS25.213, v3.2.0, 
March 2000; and 3rd Generation Partnership Project, "FDD: 
Physical Layer Procedures," 3GPP Technical Specification, 
TS25.214, v3.2.0, March 2000 (collectively "3GPF'), or in 
the CDMA2000 standard in TIA, Interim V&V Text for 
CDMA2000 Physical Layer (Revision 8.3), Mar. 16, 1999 
("TIA"), an option is to employ smart antenna technology in 
the base station. A smart antenna can suppress interfering 
signals from different direction of arrival angles (DOAs) 
from the desired users by using spatial diversity. Smart 
antenna technologies attract much attention these days as 
they support more users with a high quality of service and 
high data rates, up to 1.92 mega bits per second (Mbps). 
Examples of smart antennas and system architectures 
employing smart antennas may be found in commonly- 
owned co-pending U.S. application Ser. No. 09/610,470, 
filed Jul. 5, 2000, entitled "Smart Antenna with Adaptive 
Convergence Parameter;" Ser. No. 09/661,155, filed Sep. 13, 
2000, entitled "Smart Antenna with No Phase Calibration 
for CDMA Reverse Link;" and Ser. No. 09/669,633, filed 
Sep. 26, 2000, entitled "New Cellular Architecture." The 
contents of all of these applications are hereby incorporated 
by reference herein. 

Despite the interest in smart antenna technology 
generally, little attention has been paid to pseudonoise (PN) 
code acquisition in CDMA systems that employ a smart 
antenna at a base station. As used herein, PN code acquisi- 
tion refers to a portion of a process referred to in the art as 
synchronization. Synchronization is generally regarded as 
encompassing two processes: PNcode acquisition (in which 
a phase error for a known PN code is resolved to within a 
specified boundary — that is, a coarse PN phase code error 
correction), and PN code tracking, in which fine differences 
in PN phase code errors are detected and corrected. Thus, 
PN code acquisition, despite its misleading moniker, refers 
to a coarse correction for a PN phase error between a 
receiver and a transmitter (e.g., a base station and a mobile 
unit, or vice-versa), and does not refer to a process by which 
a PN code itself (as opposed to a PN code phase error) is 
detected. 

PN code acquisition may be difficult when the smart 
antenna weight vector does not correspond to the desired 
signal's DOA (because the smart antenna will suppress 
signals from other DOAs). Because of this potential 
problem, existing systems use only one (omnidirectional) 
element output out of the M array elements for PN acqui- 
sition purposes. See F. Adachi, M. Sawahashi, and H. Suda, 
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"Wideband DS-CDMA for Next-Generation Mobile Com- 
munications Systems/' IEEE Communications Magazine, 
pp. 56-69, September 1998 ("Adachi et al"). Thus, the 
benefit of a smart antenna (e.g., the potential gain of the 

5 smart antenna) has not been used for PN code acquisition. 
This causes another problem as the mobile units in the 
system may transmit signals with low power because of the 
expected high smart antenna gain at a base station, with the 
result that the received signal-to-interference-plus-noise- 

10 ratio (SINR) at a base station may not be sufficient for PN 
code acquisition when only one element is employed. 

BRIEF SUMMARY OF THE INVENTION 

The present invention is an efficient PN code acquisition 

15 scheme which employs multiple elements of an antenna 
array and an adaptive threshold. The invention is particularly 
useful for CDMA wireless communications, especially for 
Direct Sequence (DS-)CDMA wireless communications. 
The basic structure of preferred embodiments is the combi- 

20 nation of a conventional PN correlation searcher, an adap- 
tive beamformer, and an adaptive threshold setting circuit. 
During each observation interval, which consists of multiple 
chips, the adaptive beamformer adaptively updates the 
weight vector for the smart antenna elements (as used 

25 herein, a smart antenna element refers to a single antenna, 
such as an omnidirectional antenna, in an array of antennas 
that collectively form the smart antenna) using the accumu- 
lated received signal despread with trial PN code phase error 
as input, preferably at the chip rate. The adaptive beam 

30 former may use any one of a number of algorithms for this 
purpose. In a preferred embodiment, a normalized least 
mean square algorithm is used. A spatially correlated signal 
is then formed by weighting an accumulated value of the 
signal received by each antenna in the array over the 

35 observation period with the corresponding final weight of 
smart antenna weight vector as calculated by the adaptive 
beamformer algorithm. This spatially correlated signal is 
then compared to a threshold to determine whether PN code 
acquisition has occurred. If the threshold is exceeded, PN 

40 code acquisition is declared. Otherwise, a new trial PN code 
phase error is selected and the process is repeated. 

In conventional serial search algorithms, the above- 
mentioned threshold for PN code acquisition has been fixed 
and can be calculated from a given false alarm probability P f 

45 and a given bit-energy-to-interference power spectral den- 
sity ratio Ej/I^. But in a real mobile fading environment, 
Ej/I^ often varies. In preferred embodiments, an adaptive 
threshold setting algorithm is employed to adapt to varying 
environment for efficient PN code acquisition. An adaptive 

50 threshold setting algorithm has been analyzed for a receiver 
with a single antenna element, in Kwonhue Choi, Kyung- 
whoon Cheun, and Kwanggeog Lee, "Adaptive PN code 
Acquisition Using Instantaneous Power-Scaled Detection 
Threshold Under Rayleigh Fading and Gaussian Pulsed 

55 Jamming," The 4' A CDMA International Conference, Pro- 
ceedings Vol. H pp. 162-169, Seoul, Korea, Sep. 8-11, 1999 
("Choi et al.") , the contents of which are hereby incorpo- 
rated herein by reference. The present invention develops an 
adaptive threshold setting algorithm for a receiver with 

60 multiple array elements. The adaptive threshold setting 
circuit is actually an average power estimator in preferred 
embodiments. While the adaptive beamformer updates the 
weight vector w(i) adaptively, the power estimator estimates 
the instantaneous received signal and interference power 

65 prior to PN code despreading for each observation interval 
of NT C seconds. The average estimated power is then 
employed to scale a fixed reference threshold to create the 
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above-mentioned adaptive threshold used to determine FIG. 1 shows a block diagram of the proposed PN code 

whether or not PN code acquisition has been achieved. The acquisition scheme for a DS-CDMA system configured in 

PN code acquisition time with the proposed PN code acqui- accordance with one embodiment of the present invention, 

sition scheme with M=5 elements, by way of example, can The input signal r'(t) received from the each element 101 of 

be 210% shorter at a given SINR than the PN code acqui- s the smart antenna array is passed through a demodulator 102 

sition time for a system, such as the system described in and a matched filter 103. The output of each matched filter 

Adachi et al., in which only a single element is used for PN 103 is sampled at the chip rate by a chip sampler 104, and 

code acquisition. correlated (despread) with the trial PN code phase error by 

In one aspect of the present invention, the PN code a correlator 105. The output of the correlators 105 is stored 

acquisition system may be applied to a base station, wherein io by accumulators 106. The output of the accumulators 106 is 

the antennas are in the base station. According to another used DV an update process for the smart antenna weight 

aspect of the present invention, the system is applied to a vector. 

mobile station, wherein the smart antennas are in the mobile In preferred embodiments, the weight vector of the smart 

station. antenna is updated with an LMS algorithm. In highly 

™ T ™ ^™,™.™^ T M .„,« T 15 preferred embodiments, the normalized least mean square 

BRIEF DESCRIPTION OF THE DRAWINGS ^ MS) in § Hayldn> Adapti ^ filter 

The features, objects, and advantages of the present Theory, 3rd Edition, Prentice Hall, Upper Saddle River, N.J., 

invention will become more apparent from the detailed 1996 ("Haykin"), is employed to update the weight vector, 

description set forth below when taken in conjunction with Other vector update algorithms are also possible, including 

the drawings in which like reference characters identify 20 those algorithms discussed in commonly assigned 

correspondingly throughout and wherein: co-pending U.S. application Ser. 09/610,470, filed Jul. 5, 

FIG. 1 shows a block diagram of a proposed PN code 2000, entitled "Smart Antenna with Adaptive Convergence 

acquisition scheme for a DS-CDMA system configured in Parameter." It is not necessary for the vector update to be an 

accordance with one embodiment of the present invention; LMS algorithm. During each observation interval of NT C 

FIG. 2 shows a flow chart of the proposed PN code 25 ( where T c * the chi P P eriod and N is the number of chip 

acquisition scheme for a DS-CDMA system configured in P eriods m m observation interval) seconds, an adaptive 

accordance with one embodiment of the present invention; beamformer 107 updates the weight vector w(i) adaptively, 

* u .u .l »■ i a 1 *■ ei I preferably at the chip rate. The input to adaptive beamformer 

FIG. 3 shows the theoretical and simulation false alarm K M ■ a/ , 4 + n£ » * & j j- 

. . „ ... t . t _j. 107 is the accumulator 106 output alter despreading using 

probabilities ? f versus bit -energy-to-interference power , n . § . . nNr . . 1 , , f & & 

r . , , 1 . . ~ n . J , % t - , c t 30 the trial PN code phase error as described above, 

spectral density ratio EJ\ 0 in dB for M-l, 3, and 5 elements r 

cases under a fading environment, indicating that the PN In preferred embodiments, the vector update process 

code acquisition scheme configured in accordance with one starts the observation interval begins. This is done in 

embodiment of the present invention is robust against the order to conserve processing power. The sample time index 

received signal power variations; " iss K ( where i » the chip or sample tune index) on which the 

FIG. 4 shows the theoretical and simulation detection u P date sta f is t chosen to reflect < he time T ^ ukc6 f ° r th * 

probabilities P, versus bit-energy-to-interference power ada P b ^ e algorithm to ™ nver S c ' F ° r m 

spectral density ratio E^ 0 in dB for M-l, 3, and 5 elements ^ mb °? ime ° tS ™** th *. °T K L ™ * algorithm 

under fading environment, indicating that the PN code ***** *™ .™ d 10 whl f h [ number of chips N in the 

acquisition scheme configured in accordance with one 40 ^ to .J^£/^ * 256 > N ' * ™ ™ 2 > 

embodiment of the present invention improves detection ^ hlch provides for 256-192=64 iterations of the a gonthm. 

probability P, significantly; and ™* reflect * ^ fact that the ^ onihm *n*My ^ 

» . , . ■ , j - i • approxunately 40 iterations to converge and includes a 

FIG. 5 shows the theoretical and simulation average ^£ ^ factor 

acquisition time T„_ versus bit-energy-to-interference y . vT 

power spectral density' ratio EJ\ 0 in dB for M-l, 3, and 5 45 ^ u P date P^cess ends at i-N-2 for the first observation 

elements under fading environment, indicating that the PN in *™ 1 i The ^ wei g ht vector calculated b y 

code acquisition scheme configured in accordance with one u P datc algorithm in an observation interval which can be 

embodiment of the present invention improves average regarded as a best weight vector w,^ is used to generate an 

acquisition time T^ significantly. antenna m y s P atial correlation output z with_the last 

50 accumulation output at i«N-l. The square of z is then 

DETAILED DESCRIPTION compared to an adaptive threshold (discussed in further 

In the following detailed description, many specific detail below) to deterrmrie whether PN code acquisition has 

details, such as chip rates and numbers of elements in smart been achieved. If the threshold is exceeded, the PN code has 

antenna arrays, are provided by way of example in order to been acquired and the process ends. Otherwise, a new trial 

provide a thorough understanding of the present invention. 55 PN code P hase error * selected and ^ algorithm is 

The details discussed in connection with the preferred repeated. The weight vector adaptation is performed at each 

embodiments should not be understood to limit the present successive observation interval of NT C seconds until PN 

invention. code acquisition is achieved. 

The present invention is described below in the context of The preferred LMS algorithm takes only 5M computa- 

a reverse link from a mobile to a base station in a CDMA 60 ^ons per chip period, where M is the number of the antenna 

wireless communications system. However, the invention is array elements 101. The update rate can be sufficiently small 

equally applicable to forward link PN code acquisition. The sucn as a chip rate of 1.2288 Mcps. The present invention 

•invention is also believed to be particularly applicable to employs preferably all array elements 101 of an array for a 

Direct Sequence CDMA (DS-CDMA) systems, and thus PN correlation searcher while only one element is employed 

will be discussed in the context of the same herein, but 65 ^ the existing literatures ("Adachi et al."). 

should not be understood to be limited to DS-CDMA The following discussion explains the above-described 

systems. process in greater detail. The array antenna elements 101 are 
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assumed to be substantially identical and have the same can be regarded as a modified array response vector, com- 

response to any direction and the antenna spacing is one half pared with Ao(8) in equation (1). Assume an observation 

of the carrier wavelength. The array response vector can be interval of NT C for a PN code phase shift. Energy is 

written as accumulated every chip at each branch after PN despreading 

ojfiftu-t" "* 9 e -M**-t)*» y ^ 5 105. The accumulator 106 output at iT c for the m-th branch 

can be written as 

where 0 is the DOA from the desired user signal, and T 

denotes the transpose. The received signal at the m-th / (10) 

element can be written as >«.(') = £ **(/') 

10 f«o 

r'„(f>v^a(r)fc(f-Trjc(f-Trj^^ Mt '-V»<'>c-^ w '- 1 > fa 8 

+n m , BPr (t) (2) 

where i=0, 1, . . . , N-l for the first observation interval, and 

where P is the received signal power, a(t) and <Kt) are the m=1 2 M adaptive beamformer 107 preferably 

fading ; amplitude ; md phase, respectively, b(t) .s information „ ^ d&tin wei ht vector w(i) ada dvel at ;_ N and 

data, c(t) is the PN spreading sequence, w c and y m are the ends u £=N _ 2 for ^ observalion intervaL 

angular earner frequency and phase, respectively, T c is a /rTn _ , . . .. . , ., , 

, f . . . « • qvt i i ff , / (Those of skill in the art will recognize that it is also possible 

chip interval, and x is a PN code phase offset from a ... j4 . , 4 , UT * * 

reference, which is assumed to be a random integer uni- ^bfgm updating earlier or later, even as early as the start 

formly distributed between 0 and PN sequence period L. Hie on of the observation intervaL) The same process is applied for 

n. WO in equation (2) represents an additive white Gaus- 2 ° ^ other observation intervals until PN code is acquired, 

sian noise (AWGN) bandpass interference plus noise signal. ^ normalized LMS algorithm is employed to update the 

A fading channel with a single path is assumed for simplic- vector - ^ vector adaptation is performed 

ity. The present invention can be applied for a multipath every observation interval of NT C seconds. The input to the 

channel also. A simplified search, which shifts the trial PN 25 adaptive beamformer 107 is written as 
code sequence by T c , is considered here for simplicity. The 

phase of the PN code can be shifted in a fraction of chip xCOHyitiWO • • • y^(03T. (n) 
interval in a complete search and it takes long simulation 

time. A single dwell is assumed for simplicity. The invention The last weight vector w(N-l) in an observation interval, 

can be applied for a complete and multi-dwell PN code 3Q which can be regarded as a best weight vector w^„ is used 

searches also. The PN spreading sequence c(t) is written as to generate a spatial correlation output z with the last 

accumulation output y.(N-l). The antenna array spatial 

(3) correlation output z 109 can be written as 



c(0 



:^(0fI('-»T c ) 



35 



2 = iA/V-l)£(/V-l)/(M/V) (12) 



where c(i)e{l, -1} with equal probability, n Tc is a rectan- -_Ly Y[w x (01 

gular shape pulse with unit amplitude and duration T c . The MN lzii=i 

carrier synchronization is not assumed because noncoherent MN 

PN code single dwell serial search is employed. Also, 40 = _i_V Zy) 

assume a pilot channel is used instead of pilot symbols, i.e., MN jmi 
b(t)-l for all time as the pilot channel in cdma2000. Then, 
the equivalent baseband signal r m (t) at the m-th element after 

down-conversion with carrier frequency can be written as where H stands for Hermitian, i.e., conjugate and transpose, 

raO-^O^-xr^*"^')^-^ (4 ) 45 n{N-l^A^ op!2 • ■ • ^T, (13) 

where n m (t)=n m 1 (t)+jn m e (t) is a complex AWGN with mean 
zero and variance I 0 . The output of matched filter H*(f) 103 
is sampled at chip rate and the samples r m (i) are written as 

rJO-^(0^)^^*^" 1),fa *hnjfi (5) 

where i and E c denote the chip index and chip energy, 
respectively. After correlating with the local PN code at 105, 
r m (i) becomes 

x m (04N^a(0^-ty^^0Jc(t-T) (6) 

where x is the estimated PN code phase and 



and 

50 z/"»W«(0JK™-1)W+1. (14) 

The final decision variable Z 110 is normalized and written 
as 



55 



(15) 



♦^(O-Ym-^CO-afa-l) 8 " 1 9 CO fi0 As discussed above, this final decision variable Z is then 

All antenna element outputs can be written in a vector form a threS , ho1 ?' pr fif Ma "^ " "J?? 1 ™ ' hr f J?" * S 

as calculated by an adaptive threshold setting circuit 108, in 

order to determine whether PN code acquisition has been 

i(i^v/^a(/)c(/-T)a(<i>>+fl(0]c(/-x) (8) achieved. In preferred embodiments, the adaptive threshold 

65 setting circuit 108 estimates the average power during each 

re observation interval NT C . Referring to FIG. 1 and consid- 

flfoMe* 1 ^** 0 . . . c /4w(/ >] r (9) ering |r m (i)| 2 -x m (i)| 2 , we can get the average power P as 
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for the first observation interval, where a and <(> are con- 
(16) slants. Therefore, (7) becomes 



*J,tHrm*+-*{m-1>in ^> OS /SAM, (26) 

! m N-i 5 for the first observation interval, where <t> m is constant. 

Z E l^p^COl 2 (1) When PN Code is Synchronized 



0 This case means 



J„f>/- c<Hh4hh (27) 

j~ l 10 

For a theoretical analysis, we apply perfect weight vector 

w,^ equal to the modified array vector in equation (9) to 

If Z is larger than a scaled threshold k-F, a PN code obtain an idea] performance when PN code is synchronized. 

acquisition is declared and the tracking loop is triggered for in other words, we assume that the smart antenna can track 

a single dwell serial search, where k is the fixed reference the DOA of the desired user signal perfectly as 

threshold. Otherwise, the acquisition scheme shifts the 

locally generated PN code phase (i.e., creates a new trial PN u^rafo) (28) 

code phase error), and the searcher continues until a correct durin ^ emJ iod rf an obs6rvation interval N s iSN _ 

h C °% " acquiation is ^ when pN co(Je fa synchronized . In practice / a smart 

ac ieve j. ^ antenna cannot track perfectly the DOA of the desired user. 

The probability of Z being greater than k-P is given as ^ ideaJ case fa equation (2g) prQvides ^ upper and a 

lower bound of detection probability V d and average acqui- 
(17 ^ sition time T acq , respectively. Equation (28) implies 

25 w D ^;-<-'M^^. (29) 

Substituting equations (6), (27) and (29) into (14), we 
To ease the derivation of the false alarm probability ? f and nave 
detection probability ? d later, we rewrite equation (17) using 

new decision variables U and V as r/-^a+n y ,;-(m-l)iV+i+l, (30) 

P {Z>k-P}*>P {U>V} (18) 30 WDcre n i*s~ J ^mQ) c (i- r *)> which can be easily shown to be 

a complex AWGN variable with zero mean and variance 
w jth equal to I 0 . Substituting equation (30) into equation (12), we 

get the sample mean 

35 ~z~^W c a+^ (31) 

where 



(I MN 2 . MN \ 



t/=(l-Jfc)AflVfl 2 , (19) 
V-k(MN-l)v z , (20) 

where z is the sample mean of z y . given in equation (12), and 
v z denotes the sample variance of z, given as 

» (21) 



l MN 



MN 



is a complex AWGN variable with zero mean and variance 
equal to 1 0 /MN. Substituting equations (30) and (31) into 

Under hypothesis H, when the received PN code and the 45 et l uations < 19 > md ( 20 >' we have 
local PN code are aligned with the same phase, equation (18) u^i-k)MN\^a+n\ z , (32) 

gives the detection probability P d , and under hypothesis H 0 

when PN codes are out of phase, equation (18) gives the v-^acsmk, (33) 

false alarm probability ? f as follows: SQ w here 



JXtftMW}, (22) 
P r P r {U>V\H 0 }. (23) 



^ MN 



From equations (19) and (20), we observe that both V d and 55 
V f are equal to 1 when k<0 and both equal to 0 when k>l . is the sample variance of n,. From equation (32), we can 
So we focus on the range of k between 0 and 1 in this show the decision variable U has the non-central chi-square 
invention. distribution 20 with its conditional probability density func- 

The following is the theoretical analysis. Assume a slow tion as 

Rayleigh fading channel so that the amplitude a(i) as well 60 

as the phase <|>(i) remain constant over each observation time exp{-[u + (i -k)MNEcc?)i[{\ -*)/<>]} (34) 
NT C but are independent between the observation intervals. ?<>/(" I <*) = (i -*)/„ 
Then a(i) and 4>(i) can be written as , , — , 



a(0-a,0S/£AM. (24) 
<J>(0-4>, (25) 
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where I 0 (x) is the zero-th order modified Bessel function of 
the first kind, and a is Rayleigh-distributed with its prob- 
ability density function as 



(35) S 



where oy 2 is the average fading channel power. The prob- 
ability density function (pdf) of U under synchronization is 10 



re*p(-u//J) 
Puiiu | a)p{a)da = - , u 2; 0, 



(36) 



where 



fa(l-k)(MNE c a f 2 +I 0 ). 



(37) 



From equation (33), we can show the decision variable V 
has the central chi-square distribution with 2(MN-1) degrees 
of freedom with its probability density function as below: 



PV/(V) - ^HMfj-DyMN-DriMN - 1) 



(38) 



25 



where 



(39) 



From equations (32) and (33), we know U and V are 
statistically independent because sample mean and sample 
variance are statistically independent when sampled from a 
Gaussian distribution. Then the detection probability can be 
derived to be 



(40) 



(2) When PN Code is not Synchronized 
This case means 



c(i'-t)«»c(i-t). 



(41) 



Assume the beamforming direction of smart antenna is 
different from the DOA of the desired user signal. This is 
true in general when the PN code is not synchronized. 
Assume that the weight vector can be written as 

H^r^(* r - ^-afaH^ 1 ^ ■ • ■ «T-«(<I>) (42) 

without loss of generality. Equation (40) implies 

v>>o P <„-e~ lT * n +e- j * m , 1 ^m^M. (43) 

Substituting equations (6), (41) and (43) into equation 
(14), we have 



z y »v^nc(^T)c(i-x)e / (<r" T, *' ) +e" yTl, "« / ,(Oc(^x). 



(44) 



(45) 



50 



We can show zy is a complex AWGN variable Dq with zero 
mean and variance equal to (E c a^+l 0 ). Equation (44) can be 
rewritten as 



60 
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Substituting equation (45) into equation (12), we have 



2 = 



(46) 



which is also a complex AWGN variable with zero mean and 
variance equal to (E CO /+I 0 )/(MN). Substituting equations 
(45) and (46) into equations (19) and (20), we have 



C/-(l-*)MtfK| 2 , 



(47) 
(48) 



15 



where 



MN 



20 



is the sample variance of Uq. Similarly as subsection (1), we 
can show the decision variables U and V follow central 
chi-square distribution with 2 and 2(MN-1) degrees of 
freedom, respectively, and their probability density function 
are given as below 



exp(-M/ffl 
0 ' 



30 



^- 2) exp[-v/(2(7 3 )] 

PVOW - ^MW-ll2|MW-l>r(MV - 1) ' 



where 

35 Ml-*)(E c a/ + / 0 ), 
o 2 -A(E e o/+/ 0 )/2. 



(49) 



(50) 



(51) 
(52) 



Similarly as subsection (1), we also know U and V are 
statistically independent from equations (47) and (48), due 
40 to the fact that sample mean and sample variance are 
statistically independent when sampled from a Gaussian 
distribution. Then the false alarm probability can be derived 
to be 



45 



= (1-Jfe)^- 1 . 



(53) 



We observe that the false alarm probability P^is a function 
of k, M and N only and does not depend on E c and I 0 . Thus 
the adaptive threshold algorithm can approximately achieve 
a constant false alarm rate (CFAR). 
FIG. 2 shows the corresponding flow chart of the pro- 
55 posed PN code acquisition system with the present inven- 
tion. All the vectors in FIG. 2 have M elements, where M is 
the number of antenna elements. First, we start from the 
received sequence i(i) in block 201, which is formed by the 
output of the chip samplers 104 in FIG. 1, where i is the time 
index. Next, x_(i) (the despread signal) in block 202 is formed 
by the output of the correlators 105. Next, the y.(i) in block 
203 is formed from the accumulator 106 output. The output 
of block 203 is then used in block 204 and block 208. At 
i=N, the adaptive beamformer 107 starts updating the weight 
vector w(i). Wq-(1 1 ... 1)" in block 205 is the initial 
weight vector at i=N,. In block 204, the input to the adaptive 
beamformer 107 is £ (i), which is the normalized 
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y_(i). The output of block 204 is used in updating the weight 
vector block 206. Block 206 shows the LMS algorithm for 
the adaptive beamformer 107. As shown in block 206, an 
intermediate value w(i) of the weight vector is formed at 
each update, which is then used at the subsequent update. 
The updated weight vector from block 206 is output to block 
207. At i-N-1, the final weight vector shown in block 207 
is multiplied with the final accumulator output y_(N-l) 
shown in block 208 to generate the spatial correlation output 
109 z in block 209. Block 210 calculates the magnitude 
square of z using the output of block 209 to get the decision 
variable Z 110, which is output to blocks 212. The power is 
estimated at block 211 utilizing the outputs of block 201 and 
207. In block 212, Z is compared with the threshold scaled 
by the adaptive threshold setting circuit 108 output shown in 
block 211 and then the decision as to whether to declare PN 
code acquisition is made. 

To verify the theoretical results, simulation has been done 
with the PN code acquisition scheme configured in accor- 
dance with one embodiment of the present invention. The 
simulation parameters are described as the following. 
P^O.01, the number of antenna elements M=l, 3, and 5, 
N 5 =192, N=256, bit rate=9600 bits/second, chip rate 
«1.2288xl0 6 chips/second, and penalty factor K^-256 were 
assumed. The fixed reference threshold k was chosen to 
guarantee V f to be equal to 0.01 for the adaptive threshold 
case. k«0.0179 was chosen for one element case; k=0.0108 
for three element case; and k-0.0089 for five elements case. 
The fixed reference I 0 was chosen to make Py in equation 
(59) to be equal to 0.01 at Ef/I 0 — 15 dB for the nonadaptive 
threshold case. One sample per chip was taken for simula- 
tion. Jakes' fading model was used for fading environment. 
The average fading channel power ay 2 =l, velocity v=80 
miles/hour and carrier frequency f c =900 MHz were chosen 
for fading case. Spreading factor 128 and E i yi 0 »128E c /l < , 
were used where £ b is information bit energy. The simula- 
tion parameters are summarized in Table 1. 

TABLE 1 

Simulation Parameters 



Description 




Parameter Value 


False Alarm Probability 




P f - 0.01 


Observation Interval 




N-256 


Starting Time Index of 




N, - 192 


Updating Weight Vector 






Number of Samples per Chip 




1 sample/chip 


Penalty Factor 






Bit Rate 




9600 bits/second 


Chip Rate 




1.2288 x 10 6 chips/second 


Spreading Factor 




128 


Bit-Energy-to-Interference Power 




Eb/L - 128 EJl. 


Spectral Density Ratio 






Carrier Frequency 




f c •> 900 MHz 


Mobile User Velocity 




v = 80 miles/hour 


Average Fading Channel Power 




o/-l 


Fixed Reference M » 


1 


k = 0.0179 


Threshold k l M - 


3 


k = 0.0108 


M - 


5 


k = 0.089 


Interference Power M - 


1 


I 0 = 0.9947 


Spectral Density Ratio M «> 


3 


I 0 - 0.9981 


M o 


5 


I 0 = 0.9987 
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1 k is chosen such that the PN code acquisition yield a given P f » 0.01 for 
the adaptive threshold case. When P = 1, k becomes the fixed threshold 
for the nonadaptive threshold case. That is why k is called the fixed refer- 
ence threshold. 

FIG. 3 shows the theoretical and simulation false alarm 
probabilities P f versus bit-energy-to-interference power 65 
spectral density ratio E t /l Q in dB for M-l, 3, and 5 cases 
under fading environment with a PN code acquisition 



scheme configured in accordance with one embodiment of 
the present invention. The simulation V f for M-l, 3, and 5 
cases fluctuate very slightly around the theoretical P^O.01. 
And the simulation P/s do not increase as E 6 /I 0 increases 
while P^ does for a nonadaptive threshold scheme. This 
means the present invention can approximately achieve 
CFAR and is robust against the received signal power 
variations. 

FIG. 4 shows the theoretical and simulation detection 
probabilities ? d versus bit-energy-to-interference power 
spectral density ratio E^ in dB for M-l, 3, and 5 under 
fading environment, with a PN code acquisition scheme 
configured in accordance with one embodiment of the 
present invention. At P/=0.8, M-3 and 5 case simulation 
results show 3.5 dB and 4.9 dB improvement in SINR, 
respectively, compared with M=l case. 

FIG. 5 shows the theoretical and simulation average 
acquisition time T acq versus bit-energy-to-interference 
power spectral density ratio E^ 0 in dB for M-l, 3, and 5 
under fading environment with a PN code acquisition 
scheme configured in accordance with one embodiment of 
the present invention. At T^-IO^NT,., M«=3 and 5 case can 
improve 2.7 dB and 4.2 dB in SINR, respectively, compared 
with M-l case. Therefore, when a base station employs the 
present invention, a random access user may transmit signal 
with smaller power than the single antenna case. Interfer- 
ence can be reduced and capacity of the system can be 
increased. Or at Ej/I 0 =-10 dB, the PN code acquisition time 
with the proposed PN code acquisition scheme of M=5 
elements can be significantly shortened, e.g., 210%, com- 
pared to the existing PN code acquisition scheme of a single 
element in "Adachi et al." We observe that simulation results 
are different from the theoretical results when M^2. This is 
because the smart antenna weight vector does not perfectly 
match with the array response vector of the desired user 
signal's DOA while perfect weight vector was assumed for 
the theoretical results. 

In conclusion, a PN code acquisition scheme using a 
smart antenna and an adaptive threshold setting of the 
present invention is robust against the received signal power 
variations and could improve performance 210% in PN code 
acquisition time or at least 4.5 dB in SINR to achieve the 
same detection probability when M=5 elements are used. 

Those of skill in the art will readily appreciate that various 
alternative designs and embodiments for practicing the 
invention are possible. For example, while a single dwell has 
been described above, the present invention is equally 
applicable to other PN searches, e.g., double dwell, serial 
and complete searches. Additionally, those of skill in the art 
will recognize that the various devices described herein 
(e.g., accumulators) can be implemented in hardware and/or 
software. It is therefore to be understood that within the 
scope of the appended claims, the invention may be prac- 
ticed otherwise than as specifically described herein. 

What is claimed is: 

1. A method for acquiring a pseudo-noise (PN) code 
comprising the steps of: 

accumulating despread samples obtained from a plurality 
of elements in a smart antenna array at a sampling rate, 
the samples being despread using a PN code shifted by 
a trial PN code phase error, separate accumulations 
being maintained for each element; 

updating an antenna array weight vector based on the 
accumulated despread samples; 

forming a spatially correlated signal by weighting the 
accumulated despread samples from each element with 
a corresponding weight from the updated antenna array 
weight vector and combining the weighted samples; 
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comparing the spatially correlated signal to a threshold to 
determine whether PN code acquisition has been 
achieved; and 

selecting a new trial PN code phase error and repeating 
the accumulating, updating, forming and comparing 
steps if PN code acquisition has not been achieved. 

2. The method of claim 1, wherein the sampling rate is 
equal to a chip rate. 

3. The method of claim 1, wherein the sampling rate is 
faster than a chip rate. 

4. The method of claim 1, wherein the updating step is 
performed using a least mean square algorithm. 

5. The method of claim 1, wherein the updating step is 
performed using a normalized least mean square algorithm. 

6. The method of claim 1, wherein the threshold is fixed. 

7. The method of claim 1, wherein the threshold is 
adaptive. 

8. The method of claim 7, wherein the adaptive threshold 
is determined based upon an estimated power. 

9. The method of claim 8, wherein the process of esti- 
mating the power includes the steps of accumulating 
received signals that have not been despread from each 
element to form a plurality of spread accumulated signals, 
weighting each spread accumulated signal from each ele- 



10 



20 



24. The method of claim 1, wherein an observation period 
during which samples are accumulated is constant for each 
iteration of the accumulating step. 

25. The method of claim 1, wherein an observation period 
during which samples are accumulated varies between a 
short observation period and a long observation period. 

26. The method of claim 1, wherein an observation period 
during which samples are accumulated varies between at 
least three periods. 

27. A method for calculating a threshold used to determine 
whether PN code acquisition has occurred, the method 
comprising the steps of: 

accumulating received signals that have not been 
despread from each element of a plurality of elements 
in an antenna array to form a plurality of separately 
maintained unspread accumulated signals; 

weighting each unspread accumulated signal from each 
element by a corresponding weight from an updated 
antenna array weight vector; and 

combining the weighted unspread accumulated signals. 

28. The method of claim 27, further comprising the step 
of averaging the combined weighted unspread accumulated 
signals. 

29. The method of claim 28, further comprising the step 
of repeating the accumulating, weighting, combining and 



ment by a corresponding weight from the updated antenna 25 ave raging steps during each observation period in which the 



35 



array weight vector, and combining the weighted spread 
accumulated signals. 

10. The method of claim 9, further comprising the step of 
averaging the combined weighted spread accumulated sig- 
nal. 

11. The method of claim 1, wherein the accumulating step 
is performed for a number of samples that is one less than 
the number of samples in an observation period. 

12. The method of claim 1, wherein the updating step is 
performed at a rate approximately equal to the sampling rate. 

13. The method of claim 12, wherein the updating step is 
performed for a number of samples that is two less than the 
number of samples in an observation period. 

14. The method of claim 12, wherein the updating step is 
performed for a number of samples that is less than two less 40 
than the number of samples in an observation period. 

15. The method of claim 14, wherein the number of 
samples for which the updating step is performed is based 
upon a convergence rate of an algorithm used to update the 
antenna array weight vector. 

16. The method of claim 1, wherein the accumulating, 
updating, forming, comparing and selecting steps are per- 
formed at a base station. 

17. The method of claim 1, wherein the accumulating, 
updating, forming, comparing and selecting steps are per- 
formed at a mobile unit. 

18. The method of claim 16, wherein the base station 
forms part of a code division multiple access system. 

19. The method of claim 16, wherein the base station 
forms part of a direct sequence code division multiple access 55 
system. 

20. The method of claim 17, wherein the mobile unit 
forms part of a code division multiple access system. 

21. The method of claim 17, wherein the mobile unit 
forms part of a direct sequence code division multiple access 
system. 

22. The method of claim 1, wherein the accumulating, 
updating, forming, comparing and selecting steps are per- 
formed for a reverse link. 

23. The method of claim 1, wherein the accumulating, 
updating, forming, comparing and selecting steps are per- 
formed for a forward link. 
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threshold is used. 

30. An apparatus for acquiring a pseudo-noise (PN) code 
comprising: 

a plurality of chip samplers, each chip sampler being 
connectable to receive a signal from a corresponding 
element of an antenna array, the antenna array having 
a plurality of elements; 

a plurality of despreaders, each despreader being con- 
nected to receive samples from a corresponding chip 
sampler; 

a plurality of accumulators, each accumulator being con- 
nected to receive despread samples from a correspond- 
ing despreader; 

an adaptive beamformer, the beamformer being con- 
nected to receive accumulated, despread samples from 
the accumulators, the beamformer outputting an 
updated antenna array weight vector; and 

a PN correlation searcher, the correlation searcher being 
connected to receive the updated antenna array vector 
from the beamformer and an accumulated despread 
sample from each of the accumulators, the correlation 
searcher being configured to weight accumulated 
samples from the accumulators by corresponding 
weights from the updated antenna array vector to form 
a spatially correlated signal, the correlation searcher 
further being configured to compare the spatially cor- 
related signal to a threshold to determine whether PN 
code acquisition has been achieved. 

31. The apparatus of claim 30, further comprising an 
adaptive threshold circuit connected to supply the correla- 
tion searcher with the threshold. 

32. The apparatus of claim 31, wherein the adaptive 
threshold circuit comprises an average power calculator 
circuit. 

33. The apparatus of claim 31, wherein the adaptive 
threshold circuit comprises: 

a plurality of second accumulators, each of the second 
accumulators being connected to receive a square of an 
output of a corresponding chip sampler; 

a plurality of multipliers, each multiplier being connected 
to a respective second accumulator to receive a final 
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accumulated value, each multiplier further being con- 
figured to weight the final accumulated value by a 
weight based on a corresponding weight from the 
updated antenna array weight vector; and 
an adder, the adder being connected to each of the 
multipliers to form a sum of outputs from the multi- 
pliers. 

34. The apparatus of claim 30, wherein the threshold is 
fixed. 

35. The apparatus of claim 30, wherein the chip samplers 
are configured to sample at a chip rate. 

36. The apparatus of claim 30, wherein the chip sampler 
are configured to sample at a rate faster than a chip rate. 

37. The apparatus of claim 30, wherein the adaptive 
beamformer employs a least mean square algorithm. 

38. The apparatus of claim 30, wherein the adaptive 
beamformer employs a normalized least mean square algo- 
rithm. 

39. The apparatus of claim 30, wherein the adaptive 
beamformer updates the antenna array weight vector at a 
rate approximately equal to the sample rate. 

40. The apparatus of claim 39, wherein the adaptive 
beamformer is configured to start updating the antenna array 
weight vector at a sample that is selected based on an 
estimated number of samples required for the adaptive 25 
beamformer to converge and a safety factor. 

41. The apparatus of claim 40, wherein the adaptive 
beamformer is configured to stop updating the antenna array 
weight vector at a second to last sample in an observation 
period. 30 

42. The apparatus of claim 30, further comprising the 
antenna array wherein each chip sampler is connected to the 
corresponding element of the antenna array. 

43. The apparatus of claim 42, wherein the antenna is 
located at a base station. 

44. The apparatus of claim 42, wherein the base station 
forms part of a code division multiple access system. 

45. The apparatus of claim 42, wherein the base station 
forms part of a direct sequence code division multiple access 
system. 
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46. The apparatus of claim 42, wherein the antenna is 
located at a mobile unit. 

47. The apparatus of claim 46, wherein the mobile unit 
forms part of a code division multiple access system. 

48. The apparatus of claim 46, wherein the mobile unit 
forms part of a direct sequence code division multiple access 
system. 

49. The apparatus of claim 30, wherein the correlation 
searcher is configured to perform a simple search. 

50. The apparatus of claim 30, wherein the correlation 
searcher is configured to perform a complete search. 

51. The apparatus of claim 30, wherein the correlation 
searcher is configured to perform a serial search. 

52. The apparatus of claim 30, wherein the correlation 
searcher is configured to perform a single dwell search. 

53. The apparatus of claim 30, wherein the correlation 
searcher is configured to perform a dual dwell search. 

54. The apparatus of claim 30, wherein the correlation 
searcher is configured to perform a multiple dwell search. 

55. The apparatus of claim 30, wherein a spacing between 
adjacent elements in the array is approximately equal to one 
half of a carrier wavelength. 

56. An apparatus for calculating a threshold used to 
determine whether PN code acquisition has occurred, the 
apparatus comprising: 

a plurality of second accumulators, each of the second 
accumulators being connected to receive a square of an 
output from a corresponding chip sampler, the chip 
sampler being connected to receive a signal from a 
corresponding element in an antenna array; 

a plurality of multipliers, each multiplier being connected 
to a respective accumulator to receive a final accumu- 
lated value, each multiplier further being configured to 
weight the final accumulated value by a weight based 
on a corresponding weight from the updated antenna 
array weight vector; and 

an adder, the adder being connected to each of the 
multipliers to form a sum of outputs from the multi- 
pliers. 
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